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ABSTRACT 

Imaging and spectral observations of the Puppis A supernova remnant (SNR) with the Spitzer Space 
Telescope confirm that its IR emission is dominated by the thermal continuum emission of swept-up 
interstellar dust which is collisionally heated by the X-ray emitting gas of the SNR. Line emission is 
too weak to affect the fluxes measured in broadband observations, and is poorly correlated with the 
IR or X-ray emission. Modeling of spectra from regions both in the SNR and in the associated ISM 
show that the ubiquitous polycyclic aromatic hydrocarbons (PAHs) of the ISM are destroyed within 
the SNR, along with nearly 25% of the mass of graphite and silicate dust grains. 

Subject headings: dust, extinction — infrared: ISM — ISM: individual (Puppis A) — - supernova 
remnants 


1. INTRODUCTION 

The supernova remnant (SNR) shock processing of in- 
terstellar dust plays an extremely important role in de- 
termining the dust composition, abundance, and size dis- 
tribution, which in turn have a great effect on redistribut- 
ing starlight, regulating the temperature of most of the 
diffuse gas in the ISM, and determining the rates and effi- 
ciencies of chemical reactions in molecular clouds. Com- 
pensating for these destruction processes in the diffuse 
ISM are various gas-grain and grain-grain interactions 
that grow icy and refractory organic mantles on the sur- 
viving dust grains in the denser phases of the ISM. 

Determining the efficiency of the various grain destruc- 
tion processes is of great astrophysical importance since 
they play a major role in our understanding of the sources 
and origin of interstellar dust. Knowledge of the grain 
destruction efficiency is crucial for determining the ori- 
gin of the large amount of dust observed in the early 
Universe (Dwek et al. 2007), whether it could all be 
formed by SNe (Morgan k Edmunds 2003; Maiolino et 
al. 2004), or in AGB stars (Valiante et al. 2009; Cher- 
chneff k Dwek 2010). It also plays an important role in 
understanding the dust abundance and the origin of the 
elemental depletion pattern in the local universe (Jones 
et al. 1996; Jones 2004; Tielens 1988; Dwek 1988). 

Only in the youngest SNRs (ages < 1000 yr) is the 
emission still dominated by the ejecta of the supernova. 
Emission of most SNRs is instead dominated by the 
swept up ISM. The ideal measurement of dust processing 
in shocks is to examine a population of ISM dust before 
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and after passage through a SNR shock front. While 
many SNRs exhibit post-shock IR emission (Arendt 
1989, Saken et al. 1992, Reach et al. 2006), only of 
fraction of those have clear indication of interaction with 
specifically identifiable clouds in the surrounding ISM. Of 
these, the Puppis A SNR is probably the clearest case in 
which the cloud being swept up exhibits identifiable IR 
emission. 

Puppis A is a relatively young SNR with an age of 
~ 3700 yr (Winkler et al. 1988), and lying at a dis- 
tance of ~ 2.2 kpc (Dubner k Arnal 1988; Reynoso et 
al. 1995). It is a strong X-ray source dominated by the 
shocked ISM, but there are isolated features in the SNR 
that can be identified as SN ejecta (Hwang et al. 2008; 
Katsuda et al. 2008, 2010). The synchrotron radio emis- 
sion is similar in structure to the X-ray emission, showing 
a partial shell significantly brighter on the northeast side, 
but the two are not well correlated in detail (Castelletti et 
al. 2006). Optical emission has an entirely different mor- 
phology. The brightest features are a ridge of filaments 
along the shock front on the NW side of the SNR (Baade 
k Minkowski 1954: Goudis k Meaburn 1978). On a large 
scale, these coincide with the radio and X-ray emission, 
though small scale details differ. These filaments have 
very strong [N II] emission, indicating that circumstellar 
material is being shocked. Optical emission elsewhere 
in the SNR consists of more randomly scattered knots 
and filaments. Some of the optical knots exhibit O-rich 
spectra that characterize them as SN ejecta (Winkler k 
Kirshner 1985; Sutherland & Dopita 1995). The kine- 
matics of these ejecta knots provide age and distance es- 
timates of the SNR. Distance estimates are also provided 
by the Galactic rotation distances of observed H I and 
CO clouds that are associated with the SNR (Dubner & 
Arnal 1988: Reynoso et al. 1995: Paron et al. 2008). 

AtIR wavelengths Puppis A was first observed as part 
of the IRAS full-sky survey (Neugebauer et al. 1984). 
Puppis A was bright and well resolved from the con- 
fusion of the background ISM at 25 and 60 fin l. while 
emission from the SNR was much more difficult to dis- 
tinguish at 12 and 100 pm. At arcminute resolution the 
IR emission was observed to be well correlated with the 
X-ray emission (Arendt et al. 1990). This is perhaps the 
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best such correlation of any Galactic SNR. The bright- 
est feature at IR (and X-ray) wavelengths is a compact 
feature known as the Bright Eastern Knot (BEK). The 
BEK is at an indentation in the perimeter of the SNR 
where the SNR is in contact with a neighboring molec- 
ular cloud (Dubner & Arnal 1988), known as the “East 
Cloud” . The East Cloud is distinguishable in the broad- 
band IR imaging, which is a rarity for clouds adjacent to 
SNRs. In most other cases of SNR-cloud interaction, the 
cloud is confused with all the other clouds along the same 
line of sight. Using the IRAS data, Arendt et al. (1991) 
had compared the broadband spectral energy distribu- 
tion of the BEK with that of the East Cloud to estimate 
that as much ~ 30% of the mass of swept up dust was 
being destroyed. However, with only the limited broad- 
band SED provided by the IRAS data, there were no 
direct confirmations of the composition of the pre- or 
post-shock dust, or that the IR emission was unaffected 
by any significant line emission components. 

Thus observations of Puppis A were undertaken with 
the Spitzer Space Telescope to provide a complete high 
spatial resolution map of the SNR and the clouds with 
which it is interacting at 24, 70 and 160 pm, and to ob- 
tain 5-40 p m spectra of the interstellar dust before and 
after it is swept up by the SNR. Section 2 of this pa- 
per describes these observations and the data reduction. 
The resulting images and spectra are presented in §3, 
along with models which can be used to characterize the 
changes in the dust grain size distribution induced by the 
passage though the shock and into a hot dense environ- 
ment. In §4, we discuss the implied amount of grain de- 
struction, and the relations between the gas temperature 
and density (or ionization timescale) as derived indepen- 
dently from the IR and X-ray modeling. We summarize 
our results in §5. 

2. SPITZER OBSERVATIONS AND DATA REDUCTION 

The Spitzer observations were carried out under pro- 
gram PID 50270. The observations employed all three 
instruments on Spitzer, The Multiband Imaging Pho- 
tometer for Spitzer, MIPS, (Rieke et al. 2004), The In- 
fraRed Array Camera, IRAC (Fazio et al. 2004), and 
The InfraRed Spectrograph, IRS (Houck et al. 2004). 

2.1. MIPS 

The entire SNR and the East Cloud were observed 
using MIPS. Two Astronomical Observation Requests 
(AORs: 27514624, 27515136) with fast 1° scan legs and 
148" offsets provided complete coverage at 24 and 70 gin. 
The total integration time per pixel on the sky was 31.4 s 
and 15.7 s at 24 and 70 pm respectively. Due to the fast 
scan rate and bad pixels in the 160 pm detector array, 
there were periodic gaps in the 160 pm coverage. 

The individual frames generated from each AOR. were 
processed by the Spitzer Science Center (SSC) to produce 
the Basic Calibrated Data (BCD). The SSC also ran a 
post-BCD pipeline that mosaicked the individual frames 
of the BCD into a single mosaic for each AOR. The post- 
BCD mosaics for the 24 pm data were very good. There- 
fore our processing of these data merely consisted of re- 
projecting the 2 AORs to a common grid (2.45" pixels 
in the scan orientation) and constructing a weighted av- 
erage of the AORs using the given weights. Additional 
variants of this image were then created with pixel scales 


and resolution to match the 70 and 160 pm data, and in 
celestial coordinates. 

The post-BCD mosaics for the 70 pm data were not 
suitable for use. The 70 pm BCD frames data were repro- 
cessed using the least-squares self-calibration procedure 
of Fixsen et al. (2000). As applied here, the procedure 
represents the data at each pixel i of each frame, D l , as 
D l = S a + F p + F q , where S a is the intensity of the sky 
at location a, F p is a detector offset for detector pixel p 
(constant for all frames), and F q is another offset which 
is constant for all detector pixels but varies from frame 
to frame q. The self-calibration procedure is then ap- 
plied to derive the detector calibration parameters, F p 
and F q . As we began with the individual BCD frames, 
these detector offsets represent delta-corrections that ac- 
count for transient or temporary effects that were not 
properly removed by the standard BCD pipeline. These 
offsets were subtracted from the data D\ and then these 
corrected BCD frames were remapped to obtain the self- 
calibrated sky map S a . This procedure worked well to 
remove residual responses after the periodic stim flashes, 
and to tie together the different scan legs in a consistent 
manner. However, there was still a periodic striping ap- 
parent in the images. This striping could be mitigated 
by median filtering the data in the time domain, before 
the self-calibration is performed. The effect of the width 
of the median filter was investigated. Narrower filters 
produced cleaner results but at the expense of removing 
some of the large scale emission; wider filters were not as 
effective in removing the striping. Intermediate values 
were not a satisfactory balance. Thus, the 70 pm mo- 
saic presented here was created by Fourier transforming 
a narrow-filtered mosaic and a wide-filtered mosaic and 
combining the small-scale structure (< 4.7') from the for- 
mer with the large-scale structure (> 4.7') of the latter, 
before performing the inverse Fourier transform. The re- 
sult is far cleaner than the post-BCD mosaics. However, 
there are some linear structures in the SNR that happen 
to be aligned with the scan direction and appear to be 
filtered out by this procedure. Some of these are real, as 
they have clear 24 pm counterparts, but in other cases 
the features are artifacts that are better removed. Thus, 
some care is needed in the quantitative analysis of the 
70 pm emission. 

Significant gaps and striping in the 160 pm data pre- 
sented an even greater challenge. The standard post- 
BCD mosaics provided by the SSC were strongly striped. 
The standard filtered versions, which are also provided 
by the SSC, were free of striping, but also clearly lacked 
large scale structure in the emission as compared to the 
unfiltered versions. Both mosaics contained gaps. Thus 
the first step in our processing was to interpolate across 
the gaps in the filtered post-BCD mosaics by using a 
distance weighted average of all good pixels within a 5 
pixel radius of each missing pixel (weight = 1 /distance 2 ). 
To restore the large scale structure, a 5th order poly- 
nomial background was fit to the unfiltered post-BCD 
mosaic and then added to the “filled” version of the fil- 
tered post-BCD mosaic. This was performed on the two 
AORs independently. The AORs were reprojected to a 
common grid and combined using a weighted average, as 
was done for the 24 pm data. The general result appears 
good, with the caveat that the large scale structure may 
be inaccurate. 
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2.2. IRAC 

IRAC observations were much more limited than the 
MIPS observations, and only targeted a single 5' x 5' 
field of view at the location of the BEK where the spec- 
troscopy was to be obtained. Two pointings were re- 
quired to obtain coverage at all 4 IRAC wavelengths, so 
adjacent 5' x 5' files were serendipitously imaged to the 
north of the BEK at 4.5 and 8 pm and to the south at 
3.6 and 5.8 pm. Each set of pointings used 12-second 
frame times and a large 12-point Reuleaux dither pat- 
tern to reach a nominal integration time of 124.8 s per 
pixel, and nominal sensitivities of 0.02, 0.02, 0.07, and 
0.09 MJy sr _1 at 3.6, 4.5, 5.8, and 8 pm respectively. 
Being the brightest portion of the SNR at 24 pm, this 
would be the most likely location to detect any thermal 
emission at 3.6 - 8 pm. Also the BEK is one of the rel- 
atively few places in the SNR where there are optical 
line-emitting filaments, which may have IR counterparts 
at these wavelengths. The data from the single AOR 
(27514112) were processed using the least-squares self- 
calibration (Fixsen et al. 2000). 

2.3. IRS 

The primary IRS data were collected as two spectral 
maps using both orders of the long low module (LL1 and 
LL2; AOR 27515136) and both orders of the short low 
module (SL1 and SL2; AOR 27514368). The area cov- 
ered jointly by both modules in all four spectral orders 
is ~ 3.5' x 3.5', with additional outlying coverage in dif- 
ferent directions for different orders. The LL spectral 
map completely covered the region. However, we chose 
to undersample the SL spectral map in the cross scan di- 
rection to avoid the excessive amount of observing time 
that would be required to fully cover the desired field. 
A third AOR (27514880) used staring mode to target a 
cluster of 7 positions (3 locations in the East Cloud and 4 
background fields) with the same integration times that 
were used in the spectral maps. The data from each 
AOR were reduced and spectra for selected regions were 
extracted using CUBISM (Smith et al. 2007). However, 
the reported uncertainties were clearly smaller than the 
dispersion in the spectra as a function of wavelength. 
Therefore, the CUBISM uncertainties were renormalized 
(by factors of ~ 3) such that they matched the empiri- 
cal dispersion of the spectra. CUBISM was also used to 
produce spectral line maps for emission lines of [Ne II] 
12.8pm, [Fe II] 17.9/nn, [Fe Il] + [0 IV] 26pm, [Ne III] 
15.6pm, [S III] 18.7pm. and [S III] 33.5pm. For the 
LL spectra the background was selected from the near- 
est of the background pointings outside the SNR. at 
(a. 5) = (8 h 24 m 53.7 s , — 42°58'16") (J2000). At the SL 
wavelengths, this location was clearly too bright to use 
as a background at the location of the BEK. Thus for the 
SL spectra, the average of two more distant, but darker, 
locations were used (a, 5) = (8 h 25 m 32.9 s . — 43°07'32"), 
(8 h 25 m 23.7 s . — 42°48'24"). 

No extinction corrections were applied to the data. 
The integrated H I column density toward the BEK in 
the velocity range 0-16 km s _1 is ~ 0.15 x 10 22 cm -2 
(Dubner Arnal 1988), while X-ray modeling of the 
BEK indicates column densities from this value up to 
~ 0.4 x 10 22 cm -2 . These are fairly typical values for 
a 2.2 kpc line of sight in the Galactic plane. Using the 


extinction curve of Mathis (1990) to convert the column 
densities, we find that the most extreme extinction is 
A\ < 0.13 at the 9.7 pm silicate absorption feature, with 
Aa being several times lower at other wavelengths from 
5-40 pm. 

3. RESULTS 

3.1. Images 

A color image of the Puppis A SNR at 24 and 70 pm 
and X-rays is shown in Figure 1. The separate MIPS im- 
ages of the SNR and surrounding ISM are shown in Fig- 
ure 2. Puppis A is a large, bright, and highly-structured 
source at 25 and 70 pm, with an extent of ~ 50'. There 
is a general gradient in brightness, decreasing from the 
northeast to the southwest, which coincides with increas- 
ing distance from the Galactic plane. The bulk of the 
emission consists of filamentary structures that often ap- 
pear to be composed of shorter (~ 1') arced segments, 
which give the SNR a scalloped appearance. Many of 
the filaments appear to be unresolved (< 6") in width. 
There are only a few locations where the emission is bet- 
ter characterized as unresolved knots or clusters of knots. 
The BEK is resolved into several intersecting arced fila- 
ments and diffuse emission. This detail is shown in Fig. 
3 and identified with the nomenclature used by Hwang 
et al. (2005). These structures do not appear very dif- 
ferent from the rest of the SNR except for being several 
times brighter and marking a distinct indentation in the 
periphery of the SNR. 

At 160 pm emission is not evident from the SNR, 
including the BEK. Extrapolations of IRS spectral fits 
(§3.2.2) predict 160 pm brightnesses of ~ 6 - 12 MJy/sr 
for the BEK. The mean intensity of the 160 pm cirrus 
is ~ 130 MJy/sr with an RMS variation of ~ 15 — 20 
MJy/sr. Thus, as most of Puppis A is many times fainter 
than the BEK, it is no surprise that the SNR is not vis- 
ible at 160 pm. 

Comparison with X-ray images in Figures 2 and 3 
shows a strong correlation between IR and X-ray emis- 
sion on all spatial scales down to the ~ 6" resolution 
of the 24 pm MIPS data. In comparison, the correla- 
tion with the radio emission (Castelletti et al. 2006) 
is relatively poor, especially on small spatial scales. 
The distribution of optical filaments and knots in Pup- 
pis A is very different from that of the IR, X-ray, or 
radio emission. Indeed, only some of the brightest 
optical features have detectable 24 pm counterparts. 
The most notable of these is the O-rich ejecta of the 
“Omega Filament”, which is the only optical filament 
to have a bright easily-distinguished 24 pm counterpart 
(a, 5) = (8 h 22.8 m , -42°46.6'). The adjacent X-ray emit- 
ting ejecta observed by Katsuda et al. (2008) also has 
highly correlated 24 pm emission, but it is three times 
fainter than that of the Omega Filament. 

The IRAC imaging only covered the region of the 
BEK. Emission associated with the BEK is very faint 
and confused. Therefore the IRAC images shown in Fig- 
ure 3 have had stellar sources removed by subtraction 
of scaled versions of the 4.5 pm image for the 8 pm 
data, and scaled versions of the 3.6 pm image for the 
5.8 and 4.5 pm data. This subtracts stellar emission 
which is relatively blue in color, with minimal alter- 
ation of the relatively red emission of the BEK and the 
ISM. At 3.6 pm such a template subtraction is not pos- 
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Fig. 1. A composite image of Puppis A as seen at 70 /an (red), 24 p m (green), and X-rays (blue). Most SNR features are seen in all 
three bands. There is some intrinsic variation in the relative colors, but structure in the 70 pm background also leads to some apparent 
color features. 


sible, so the version shown is generated from the point 
source fitting and subtraction when processing the BCD 
data with MOPEX/APEX (Makovoz k Marieau 2005). 
The BEK is most clearly seen at 8 pm where emission 
lies along the leading (eastward) edges of the brightest 
24 pm features (the "Compact Knot” and the "Bar”). 
The peak intensity is 1(8 pm) < 0.3 MJy/sr. The ratio 
7(8 pm)/7(5.8 pm) = 2.7. The ratios 1(8 pm)/7(4.5 pm) 
and 1(8 pm)/7(3.6 pm) appear to be ~ 10. These IRAC 
colors are consistent with PAH emission as seen at other 
locations in the ISM, e.g. in the direction of the Galactic 
center (Arendt et al. 2008). 

3 . 1 . 1 . Global Properties 


The MIPS observations cover all of Puppis A and hence 
we used them to measure the total flux density of the 
SNR at 24 and 70 pm, the main uncertainty being con- 
fusion from ISM along the line of sight. To minimize this 
problem, we used the 160 pm image, where we do not see 
the SNR, as a template for the ISM emission. This was 
scaled and subtracted from the 24 and 70 pm images. 
The resulting maps had a smoother flatter background 
with a mean value near 0. The SNR intensity was in- 
tegrated within an irregularly shaped aperture that was 
defined jointly by surface brightness thresholds of the 
24 and 70 pm images. An expanded annular aperture 
around this region was used to fine tune the local back- 
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Fig. 2. Spitzer MIPS images of the Puppis A SNR and the surrounding ISM at 24, 70 and 160 /im. The ROSAT HRI image is shown 
in in the lower right panel. On the 24 pm image, contours of CO emission (12-17 km s _1 ) of the East Cloud (Dubner & Arnal 1988) are 
superimposed, and a circle marks the Bright Eastern Knot (BEK) where the SNR interacts with outlying clumps of this cloud. At 24 and 
70 pm, IR emission of the SNR is well-correlated with the X-ray emission, but at 160 pm it is lost in the confusion of the ISM. 


ground level at the SNR. The integrated flux changed by 
~ 5% depending on the width of the annular aperture 
used to set the local background. This was used as the 
quoted uncertainty of the integrated flux densities. Two 
upper limits for the 160 pm flux density were measured. 
The first simply integrated all the emission within the 
area of the SNR after subtraction of a flat background 
at the local minimum intensity (90 MJy sr _1 ). This 
limit includes much of the ISM along the line of sight. 
A more restrictive limit was obtained by assuming that 
the 160 pm emission is spatially correlated with the 24 
and 70 pm emission. The slopes of the linear correlations 
between the 160 pm intensity and the 24 and 70 pm in- 
tensities were used to scale the shorter wavelength flux 
densities to a 160 pm flux density. The resulting flux 
densities are listed in Table 1. The results are very sim- 
ilar to the IRAS measurement of Arendt et al. (1991). 


but with smaller uncertainties resulting from the better 
spatial resolution and the use of the 160 pm data to re- 
duce uncertainties associated with the line of sight ISM. 

3.2. Spectra and Models 
3.2.1. The East Cloud 

The Spitzer IRS spectra of the East Cloud are impor- 
tant, to this study in that they allow us to characterize the 
composition and grain size distribution of the interstellar 
dust within a cloud that is physically associated with the 
SNR. This provides a direct measure of the state of the 
dust prior to its processing by the SNR shock. For most 
other SNRs. the nature of the interstellar dust must be 
assumed or simply measured only as a line of sight inte- 
gral. 

The primary location for the East Cloud spectrum was 
selected near the brightest portion of the cloud, near 
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Fig. 3. The top panels show the full resolution 24 pm image and the Chandra ACIS image at the location of the BEK. The other 
panels show the corresponding Spitzer IRAC images. The stellar sources are subtracted to emphasize the faint diffuse emission. (See text 
for details.) The Compact Knot and Bar are both detected, but are faint and confused at 3.6 - 8 pm. 
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TABLE 1 

Integrated Flux Density of Puppis A 


Wavelength 

(pm) 

Flux Density® 

(Jy) 

Satellite 

12 

< 100 

IRAS 

24 

230 ± 11 

Spitzer 

25 

245 ± 120 

IRAS 

60 

1195 ±350 

IRAS 

70 

1025 ± 51 

Spitzer 

100 

1265 ± 990 

IRAS 

160 

< 330 b 

Spitzer 

160 

< 6000 c 

Spitzer 


a IRAS results are from Arendt et al. 

(1991), Spitzer results are from the present 
work 

b Limit for emission that is spatially cor- 
related with 24 and 70 pm emission. 

0 Limit for all emission across the extent 
of the SNR. 

(a, S) = (8 h 26 m 50.50 s , -42°55'26.0") as indicated by the 
IRAS 25 pm images. (The slit locations are indicated 
in Figure 4.) The SL and LL spectral extractions are 
necessarily over different regions due to the slit sizes and 
orientations. To maximize the signal to noise ratio, we 
extracted the average spectrum over the entire subslit 
for each spectral order. The LL and SL spectral orders 
are normalized such that the regions of spectral overlap 
coincide in intensity. The resulting spectrum is shown in 
Figure 5a. At these wavelengths, the East Cloud spec- 
trum is dominated by the typical PAH emission bands. 
The H 2 0-0 S(0) and S(l) rotational lines are detected 
at 28.2 and 17.0 /im. Thermal emission from larger dust 
grains is of increasing significance at A > 20 pm. 

We fit these spectra using a Zubko et al. (2004; here- 
after ZDA) BARE-GR-FG model, radiatively heated by 
an ambient interstellar radiation field (ISRF). The ZDA 
models are a set of models for the dust grain size (or 
mass) distributions that are derived to simultaneously 
fit average interstellar extinction, emission, and abun- 
dances. Different models are distinguished by different 
choices of grain compositions and different sets of abun- 
dance constraints. The BARE-GR-FG model comprises 
PAHs, silicate and graphite grains without any mantles, 
and is fit to elemental abundances determined from F and 
G stars (Sofia & Meyer 2001). The BARE-GR-FG model 
is the best of the ZDA bare gain models. The primary 
free parameter of our fit to the East Cloud spectrum is 
the strength of the ISRF. \ = I East Cloud) hocaU which 
was varied over the range log(\) = [—2, +3] in steps of 
0.5, with Ijocai defined by the local ISRF as characterized 
by Mathis et al. (1983). The overall normalization of 
the spectrum (proportional to the column density of the 
emitting material) is the second free parameter allowed 
in the fitting. The nominal best fit is found for a moder- 
ately enhanced ISRF. log (yj == 1. and is shown in Figure 
5a. The enhancement could be related to emission from 
the SNR itself, similar to the precursor observed outside 
of SNR N132D (Morse et al. 1996). However, the spec- 
trum is at a location that is ~ 20 pc (in projection) from 
nearest edge of the SNR. and there is no indication of a 
gradient in the ISRF with respect to distance from the 
SNR. 

The errors in the fit shown in Figure 5a, appear to be 
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Fig. 4. (top) Wide view' of the IRS coverage showing the 

BEK and the two pointings at the brightest part of the East 
Cloud and a northern extension of the cloud. A single con- 
tour at 20.5 MJy sr -1 helps delineate the extent of the East 
Cloud, (middle) IRS coverage of the Bright Eastern Knot 
(BEK) superimposed on the 24 pm MIPS image. The solid 
red and yellow lines indicate the regions of the fully sam- 
pled LL1 (19.5 — 38.0 pm) and LL2 (14.0 — 21.3 pm) spectral 
maps, respectively. The dashed red and yellow' lines indicate 
the regions of the undersampled SL1 (7.4—14.5 pm) and SL2 
(5.2 — 7.7 pm) spectral maps, respectively. The white boxes 
indicate the regions where we have extracted mean spectra for 
detail analysis. From northeast to southwest these are desig- 
nated: “Compact Knot", “Bar-N", “Bar-S”. (bottom) The 
same coverage and analysis regions are superimposed on an 
optical emission line image: Ha = red, [S II] = green, [O III] 
= blue. Only the Compact Knot region contains strong line 
emission. 
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Fig. 5. — (a) ZDA model (red line) fit to the spectrum of the 
East Cloud, (b) Spectrum of the East Cloud fit to the ZDA model, 
allowing independent normalization of each spectral order (labelled 
along the bottom of the plot), (c) Spectrum of the East Cloud 
at the second (northern) location fit to the ZDA model, allowing 
independent normalization of each spectral order. 

systematic within the different spectral orders in which 
the data are obtained. This suggests that the data in the 
different spectral orders may not be accurately normal- 
ized. Therefore, we applied an alternate strategy of nor- 
malizing each spectral order of the data independently 
to optimally match the model(s). The switch from 1 
to 4 normalization parameters does not provide any ad- 
ditional physical significance, but it does permit com- 
pensation for systematic errors between the data in dif- 
ferent spectral orders. This strategy again finds a best 
fit for log (x) = 1 as shown in Figure -5b. In this case 
the normalization used for the spectral orders [SL1, SL1, 
LL2, LLlj are [0.63,0.72,1 .35,0.82] times those used to 
normalize the orders by their overlap. As a check on 
this procedure, we similarly fit the spectrum at a second 
location, in a northern extension of the East Cloud at 
(a, <5) = (8 h 26 m 04.60 s , — 42°44 , 01.6"). These data, while 
slightly fainter, are a slightly better match to the ZDA 
model (Figure 5c), again with log (a) = 1.0. 

These results show that the ZDA B ARE-GR-FG model 


provides a good representation of the dust giving rise to 
the emission of the East Cloud. The results indicate a 
moderate (~ lOx) enhancement in the strength of the 
ISRF in the region. This is largely a result of the rising 
intensity at wavelengths > 20 /xm. As the ISRF strength 
increases, large dust grains are driven to higher temper- 
atures and the slope of this rise increases. The most sig- 
nificant discrepancy between the data and the model is in 
the region of the LL2 spectral order, 14 - 21 /xm. In this 
spectral region, it appears that the ZDA model’s absorp- 
tion coefficients (most likely for the PAHs) are roughly a 
factor of 2 too large (as implied by Figure 5a). An alter- 
nate possibility is that the model is a very good match, 
and there are systematic errors in the data that cause 
apparent discontinuities at the ends of the SL1, LL2, 
and/or LL1 spectral orders (as in Figure 5b, c). 

3.2.2. The Bright Eastern Knot 

Having demonstrated that the ZDA dust model pro- 
vides a reasonable fit to the East Cloud, we now apply 
the same model to the spectra of the three brightest por- 
tions of the BEK. The selected regions are indicated in 
Figure 4. The easternmost region corresponds to the 
“Compact Knot” feature of Hwang et al. (2005), while 
the other two cover the “Bar” region. We designate these 
latter two regions as “Bar-N” and “Bar-S” . The spectra 
are shown in Figure 6. These spectra are very different 
from those of the East Cloud: the strong PAH features 
are absent, and the peak of the thermal emission of the 
large grains is shifted from ~ 100 /xm to ~ 50 /xm. 

For the spectra shown in Figure 6, atomic line emission 
is only evident in the compact knot region. Examination 
of the spectral data cubes and emission line maps (Fig. 
7) reveals that the line emission is not correlated with 
the dust. Instead the lines are correlated with optical 
knots and filaments scattered in the vicinity of the BEK. 
Lines of singly ionized species ([Ne II] 12.8 /xm and [Fe II] 
17.9 and 26.0 /xm) are correlated with optical Ha and 
[S II] emission, while more highly ionized species ([O IV] 
26.0 /xm, [Ne III] 15.6 /xm, and [S III] 18.7 and 33.5 /xm) 
are correlated with optical [O III] emission. (We note 
that part of the optical line emission, especially the long 
linear [O III] structures in Figure 7, are part of the Vela 
SNR which lies in the foreground of Puppis A. Despite 
the superposition of the two SNRs, we see no evidence 
for any IR line or continuum emission from the Vela SNR 
that would affect the analysis of Puppis A.) No molecu- 
lar H 2 emission lines are observed in these background- 
subtracted spectra, although line-of-sight H 2 emission is 
fairly uniform over the entire spectral map prior to the 
background subtraction. This lack of H 2 is consistent 
with results from Paron et al. (2008) who find no 12 CO 
emission directly associated with the BEK. despite its 
presence in more distant portions of the East Cloud and 
even in small clumps that appear unshocked and in the 
foreground to the BEK. They conclude that either the 
gas currently being swept up is not part of a molecular 
clump, or the molecular gas is being dissociated by a 
radiative precursor. 

The interstellar dust that has been swept up into the 
SNR resides in a vastly different environment than the 
general ISM. The strong correlation of the 24 and 70 /xm 
emission with the X-ray emission indicates that the dust 
is located within the shocked X-ray emitting gas. In 
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Fig. 6. — Spectra of the BEK regions fit to the ZDA model. 
The solid red lines show the best fitting models. The dotted red 
lines show models for the same gas parameters ( T e , n e ) but with no 
sputtering (A a = 0). The circles indicate the MIPS measurements 
at 24 and 70 Jim. These are normalized to match the IRS data at 
24 /an. The 70 pm data are well fit by the models even though the 
models were fit only to the IRS data. 

such an environment, collisions with electrons are the 
dominant heating source for the dust, even if the local 
radiation field is elevated by factors of 10 as indicated 
by the East Cloud. Thus, instead of the radiation field 
strength, the free parameters of the dust model are the 
electron density, n e , and temperature, T e , of the X-ray 
emitting gas. A third free parameter characterizes the 
amount of grain destruction (if any) that has taken place 
via thermal sputtering of the dust population. Indepen- 
dent of the grain radius, a. sputtering (mainly by pro- 
tons) acts to reduce the size of the grains by an amount. 
Aa ~ n f{t ~ n e t where t is the amount of time the dust 
has been subjected to the sputtering. 

As with the East Cloud, we calculated a grid of mod- 
els, now involving a 3-dimensional parameter space. 
{n e ,T e ,Aa}, instead of only one dimension {\}. The 
gas parameters are varied logarithmically over the range 
of log 2 [n e (cm“ 3 )] = [0, 5] in ten steps of 0.5, and 
log[T e (A')] = [5.75, 7.251 in seven steps of 0.25: i.e. 
1 < n e < 32 cm“ 3 and 0.56 x 10 6 < T e < 17.8 x 10 6 


K. We varied the sputtering linearly over the range of 
Aa = [0, 0.01] pm in steps of 0.0005 pm = 5A, We apply 
the same Aa for all grain compositions. This results in 
a grid of 10 x 7 x 21 models. 

Again, each spectral order was normalized to fit each 
model in the grid to determine the single best-fitting 
(minimum y 2 ) model. However, the SL2 data (5-7.7 /um) 
were not included due to a lack of any detectable sig- 
nal. Additionally, the SL1 data (7.4-14.5 pm) showed 
evidence of an over-subtracted background: there was a 
strong negative feature at the location of the 11.3 pm 
PAH band, and the spectra show a steady slope through 
zero and into negative intensities. This cannot be easily 
alleviated by selecting a different background region be- 
cause the large scale gradient in ISM intensity (seen best 
at 160 pm, Fig. 2) and small-scale local variations hap- 
pen to place the BEK in a low intensity region relative 
to the available off-source background pointings (better 
illustrated at 24 and 70 pm, Fig. 2). Therefore, rather 
than fitting with just a single normalization coefficient, 
I data = AI mo dei, the fit involved determining the coeffi- 
cients A, B, and C in I da ta = AI mode i +BI IS m + C. The 
new B coefficient directly accounts for over-subtraction 
of the ISM background ( I ism is given by the model fit to 
the East Cloud), and the new C term allows for a more 
generic, constant error in the background. 

The fitted spectra for the three BEK regions are shown 
in Figure 6. In the figure we have also overplotted the 
mean intensities for these regions as extracted from the 
MIPS 24 and 70 pm images (open circles). These MIPS 
intensities were normalized to match the IRS data (not 
the model) at 24 pm, but were not used as a constraint 
in fitting the models. The corresponding n e , T e , and 
Aa parameters are listed in Table 2, along with the re- 
duced xt values of the fits. For comparison, the mean 
gas parameters derived from modeling the X-ray spectra 
(Hwang et al. 2005) are also listed. The IRS spectra of all 
three regions are dominated by continuum emission from 
dust that is much warmer than the typical ISM. Larger 
0.1 pm grains are at equilibrium temperatures of T ~ 48 
K, while smaller gains are radiating at higher (and fluc- 
tuating) temperatures. Thus the resulting spectra are 
warmer, and not well-characterized by a single tempera- 
ture. PAH emission is not evident in the spectra, despite 
the apparent presence of weak PAH emission in the IRAC 
imaging. The IRAC intensities indicate that the PAH 
emission is < 0.3 MJy sr -1 , which would be lost in the 
noise and systematic errors of background subtraction in 
the IRS spectra. Figure 6 also shows model spectra for 
the best fit n e and T e , but with the requirement that 
no grains were destroyed, i.e. Aa = 0.0. These model 
spectra (dotted lines) are based on the same dust content 
as the ISM, but the collisional heating within the SNR 
raises the grain temperatures which reduces the relative 
strength of the PAH features (cf. Figure 5). The exces- 
sive emission at short wavelengths is a result of excluding 
dust sputtering in these models. 

Figure 8 shows 2-D slices through the 3-D y 2 param- 
eter space at the location of minimum \- 2 . Models with 
xl < 2 min (y 2 ) are reasonable fits to the data. The 
grid of models used was fairly coarse, but the broadness 
of the xl contours in these slices indicate that higher 
precision is not demanded. The slices demonstrate two 
fundamental correlations of the parameters of the fits. 




10 


ARENDT ET AL. 


C 

.2 

O 

c 

13 

0) 

03 


c 

.2 

o 

.2 

__ 

<D 

Q 


c 

.2 

o 

.2 

o 

Q 


- 42 ° 


-43° 


-42° 


-43° 


-42° 


-43° 



Right Ascension 


Right Ascension 


Right Ascension 


Fig. 7 .— IR emission line images are compared to optical emission line images in the BEK region. IR lines from singly ionized species are 
well correlated with optical Ha and [S II], Lines from doubly ionized species correlate well with optical [O III] emission. The 26 p, m line 
likely contains both [O IV] and [Fe II] because it correlates with all three of the optical lines. The [Ne II] spectral map is undersampled by 
a factor of 2 in the SE-NW dimension, and linear interpolation has been used to estimate the missing data. 


The first is that the heating rate of a dust grain is pro- 

o /o 

portional to n e T e when the electrons are completely 
stopped by the grains (Dwek et al. 2008). This depen- 
dence is the product of the collision rate ~ n e v e [where 
v e ~ (kTe) 1 / 2 } and the kinetic energy of the electrons 
~ kT e . This is the major trend in the left column panels 
in Figure 8 and is shown as a dotted line. The spectrum 
of a large grain at its equilibrium temperature will de- 
pend on n e T e '“, and n e and T e cannot be determined 
independently. However, in practice, this degeneracy is 
broken and a minimum in \Z is localized because the 
models (and the SNR) include a range of grain sizes, 
and only the largest grains reside at equilibrium tem- 
peratures determined strictly by n e T 3/,2 . Smaller grains 
have fluctuating temperatures, and may not completely 
stop the electrons. Because of these effects, the spectrum 
is sensitive to the grain size distribution as well as the 


gas temperature and density. The second correlation is 
the inverse relationship between the electron density (n e ) 
and the amount of sputtering (Aa), as illustrated in the 
right column panels of Figure 8. Increases in electron 
density raise the heating rate and thus the grain temper- 
ature, which produces more short wavelength emission. 
Conversely, increasing the amount of sputtering reduces 
the abundance of the smallest grains in the size distribu- 
tion, leading to a decrease in the short wavelength emis- 
sion because the small grains are hotter than the larger 
grains. The center column panels of Figure 8 exhibit a 
superposition of these two fundamental trends. 

3 . 2 . 3 . The Total SNR 

As a demonstration of the usefulness of the spec- 
troscopic data versus spectral energy distributions con- 
structed from broad-band imaging, we have also fit the 
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TABLE 2 

IR Spectral Model Parameters for BEK Regions 


Region 

Aa (pm) 

n e (cm 3 ) 

T e (K) 

xl 

X-ray 

n e t (cm~ J s) T e (K) 

Compact Knot 

0.0015 

4.0 

3.2 x I0 6 

2.79 

1.70 x 10 11 

5.8 x 10® 

Bar-N 

0.0020 

5.7 

3.2 x 10 6 

2.88 

2.95 x 10 11 

7.1 x 10® 

Bar-S 

0.0025 

5.7 

3.2 x I0 6 

9.56 

2.77 x 10 11 

7.4 x 10® 



1 10 1 0 6 1 0 7 1 10 
n o (cm -3 ) T e (K) n„ (cm -3 ) 


Fig. 8.— Slices of x 2 for fitting the Compact Knot, Bar-N, and 
Bar-S spectra (top to bottom) using the ZDA dust model. Each 
2-dimensional slice passes through the minimum value of x 2 in 
the 3-dimensional (T e , n e , A a) parameter space. Contour levels 
= [1.2, 2, 4,8, 16,32] * min(x 2 ). For the Compact Knot and Bar-N 
regions, the black circles identify the projection of the minimum 
X 2 onto these slices, when the T e and Aa are constrained by the 
results of the X-ray modeling of Hwang et al. (2005). 

spectral energy distribution of the entire SNR (Table 1) 
using the grid of collisionally heated dust models. Fig- 
ure 9 shows the results. The fits are almost entirely con- 
strained by the 24 and 70 pm MIPS data points, with 
their smaller uncertainties. There is a wide range of mod- 
els that can fit the observed 24 to 70 pm color and the 
IRAS data points as well. Acceptable models can be 
found up to the tested limits (and presumably beyond) 
of the {T e .n e , Aa} parameter space. The lack of use- 
ful measurements at A < 24 pm or any constraint on 
the slope of the spectrum at 24 pm makes it impossi- 
ble to assess the relative abundance of small grains and 
the amount of sputtering. If we assume that there is no 
separate cooler component of the dust that is not cor- 
related with the 24 and 70 pm emission (i.e. adopt the 
lower of the 160 pm upper limits), then the integrated IR 
flux of the SNR is 9.0 x 10 - 8 erg s" 1 cm -2 , correspond- 



Fig. 9. Integrated flux densities of the entire Puppis A SNR. 
Solid data points are Spitzer MIPS observations. Open data points 
are IRAS data from Arendt et al. (1991). The higher upper limit 
at 160 includes all emission across the SNR, most of it likely to 
be ISM. The fainter 160 pm upper limit applies only to emission 
that is correlated with the 24 (and 70) pm emission. The red line 
indicates the best fitting ZDA dust model, while the blue lines 
indicate all the other models from the grid in {T e ,n e ,Aa} which 
have a reduced x 2 < 1. 

ing to a luminosity of 1.4 x 10 4 Lq and a dust mass of 
0.25 M© at a distance of 2.2 kpc. The luminosity and 
mass are in very good agreement with those previously 
reported (Arendt et al. 1991), although the flux is 1.32 
times higher than reported by Dwek et al. (1987). Ac- 
cordingly, the ratio of IR to X-ray emission of Puppis 
A is also higher than that of Dwek et al. (1987), with 
IRX — F//j/Fo. 2 - 4 .okeV = 6.9, which is still a factor of 
~ 5 below the theoretically expected ratio. This factor is 
too large to attribute to the grain destruction by sputter- 
ing and must represent a relatively low dust-to-gas mass 
ratio in material being swept up as concluded by Dwek 
et al. (1987). 

4. DISCUSSION 

4.1. Grain Destruction in Shocks 

The excellent correlation between IR and X-ray emis- 
sion revealed by the high-resolution MIPS imaging 
strengthens the direct association of the dust with the 
X-ray emitting gas, which had been indicated by the 
IRAS data (Arendt et al. 1991). Puppis A had not been 
observed with ISO and falls below the sensitivity limits 
of the MSX survey. The IRS spectra confirm that the 
broadband emission is dominated by thermal continuum 
emission of dust rather than atomic or molecular emis- 
sion lines. The steep drop in the continuum at shorter 
wavelengths is a clear indication that small dust grains, 
which are directly confirmed to be present in the sur- 
rounding ISM, are absent within the shocked gas of the 
BEK. The comparison of the dust grains size distribu- 
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Fig. 10. - Mass distribution of ZDA dust grains, before (black 
lines) and after sputtering by Aa = 0.002 fim (red lines). The un- 
sputtered mass distribution is that used to model the East Cloud. 
The sputtered mass distribution is appropriate for the Bar-N re- 
gion. This amount of sputtering destroys 100% of the PAHs, and 
removes 24% and 23% of the total mass in graphite and silicate 
grains respectively. 



Aa (yum) 

Fig. 11. — Fractional mass of ZDA dust grains (BARE-GR-FG 
model) eroded as a function of the total amount of sputtering Aa. 

tion [/(a) = dn/da] in the ISM with the post-shock size 
distribution shows that ~ 25% of the total dust mass has 
been destroyed in the BEK. This is depicted in Figure 
10 where the pre-shock and post-shock distribution of 
a 4 /(a) (oc a dm/ da. with dm/ da = the mass distribution 
of the dust grains) is plotted as function grain size, a, for 
each of the physical components (PAH, graphite, silicate) 
of the ZDA BARE-GR-FG model. The plot provides a 
good illustration of the fact that most of the mass of the 
dust resides in the largest dust grains despite the greater 
numbers of small grains. It also shows that under the 
action of sputtering, most of the mass loss occurs in the 
destruction of the smallest dust grains rather than the 
erosion of the largest grains. This is also shown in Fig- 
ure 1 1 which demonstrates that the dust mass destroyed 
jumps quickly as sputtering begins and then increases 
more slowly as sputtering continues to larger values of 
A a. 

Qualitatively similar results are obtained in a non- 
radiative shock examined by Sankrit et al. (2010) in 


the Cygnus Loop. There, with slightly lower post-shock 
electron densities and substantially lower temperatures, 
roughly 35% of the dust mass is destroyed over a spa- 
tial scale of ~ 0.14 pc. This spatial scale is not well- 
resolved in Puppis A, which is 4 times more distant than 
the Cygnus Loop. 

In principle, the destruction of the dust may lead to 
an increase in the strength of emission lines associated 
with atoms that are highly depleted in the gas phase of 
the general ISM, e.g. Si, and Fe. Emission of [Fe II] 
is detected at 17.9 and 26.0 /i,m, however this emission 
is spatially correlated with optical filaments and knots 
visible in Ha and [S II] rather than the dust and X-ray 
emission. There are such structures near the Compact 
Knot, but none are present near the bright parts of the 
Bar. No emission from the [Si II] 34.8 gm line is ob- 
served. This line is ubiquitous in the raw spectra, but 
then disappears completely when the background sub- 
traction is applied to the data. The lack of IR emission 
lines from the sputtered atoms is likely due to the fact 
the atoms are quickly ionized beyond the singly ionized 
stages which have strong lines in the range of the IRS 
spectral coverage. Itoh (1989) has modeled the time- 
dependent ionization of atoms sputtered from dust. His 
model 65D matches the T e = 3 x 10 6 K environment of 
the BEK. This model indicates that most of the sputtered 
Si atoms proceed beyond single ionization after n e t > 1 
cm - ' 3 yr. Thus, at the density of the BEK, freshly sput- 
tered Si atoms can only produce the [Si II] 34.8 /im line 
for < 0.25 yr. 

The spectra of the Compact Knot and the Bar re- 
veal no PAH emission, even after we implemented a fit- 
ting procedure that can recover over-subtraction of the 
background. Yet the IRAC imaging suggests weak 3.6- 
8 /rm emission at the eastern (leading) edges of the bright 
24 /rm and X-ray emission. The broadband colors of this 
emission are similar to those produced by PAHs, and 
the low surface brightness is consistent with the non- 
detection of PAHs in the spectra. The morphology and 
relative faintness compared to the 24 /im emission sug- 
gest that the IRAC data are tracing the zone behind the 
shock front (or the precursor in front of the shock) in 
which the PAHs are being destroyed. This is consistent 
with the fact that PAHs should be completely destroyed 
in v > 120 km s” 1 shocks (Micelotta et al. 2009). 

4.2. Characterization of the X-ray Emitting Gas 

In modeling the IR emission of collisionally heated 
dust, we have chosen to treat the gas properties as free 
parameters, rather than constraining them with results 
from modeling the X-ray spectra (e.g. Hwang et al. 
2005). Thus we obtain an independent characterization 
of the hot post-shock gas. The models of the X-ray spec- 
tra are characterized by the gas temperature (T e ), an 
“ionization timescale” (n e t) which indicates the length 
of time since the gas was shocked via its approach to 
ionization equilibrium, and a line of sight column den- 
sity (Nh, not an intrinsic property of the SNR). The gas 
temperature is a direct parameter of the IR modeling, 
and the ionization timescale is proportional to A a via 

n e t rifjt « 6 x 10 12 Aa (1) 

where the proton and electron densities nn, n e , are in 
cm” 3 , t is in s, and A a is in /mi, for dust that is thermally 
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and kinetically sputtered in hot gas (T e > 10 6 K) behind 
fast shocks (v s > 500 km s” 1 ) (e.g. Dwek et al. 2008). In 
detail, the coefficient of this relation is dependent on the 
grain composition (i.e. sputtering yield), and whether 
or not kinetic sputtering is acting in addition to thermal 
sputtering. 

The 3.2 x 10 6 K gas temperature that we derive in all 
three portions of the BEK is somewhat lower than cor- 
responding X-ray-derived values of T e « 5 — 7.4 x 10 6 
K across the Compact Knot and Bar regions (See Table 
2). The inferred values of Aa are significantly lower than 
the values that are indicated by the ionization timescale, 
n e t , of the X-ray modeling. Using Eq. (1) with the IR- 
derived sputtering and gas density allows us to estimate 
the time since the dust was swept up as t ss 70 — 90 
yr. This age is far shorter than the ~ 2000 yr age de- 
rived by Hwang et al. (2005) based on comparison of 
the X-ray spectrum and morphology to shocked clouds 
models calculated by Klein et al. (2003). The IR sput- 
tering timescale is also shorter than the ~ 600 — 1200 
yr timescale derived by dividing the X-ray-derived ion- 
ization timescale by the IR-derived gas density. Here we 
assume that the mean n e t is half of the maximum n e t 
derived from the X-ray spectral models, which are for 
plane-parallel single-temperature shocks and hence inte- 
grate the emission over a range of n e t. 

The remaining discrepancy suggests that the dust sput- 
tering rate of Eq. (1) is overestimated, either because 
kinetic sputtering plays a smaller than expected role, or 
because the sputtering yields are smaller. Additionally, 
the ages derived from dust sputtering may differ from the 
X-ray-derived ages because the interaction with the East 
Cloud is an ongoing process, not an instantaneous event 
that occurred at a specific date. In detail, the models 
should be integrated over time, weighted by the rate at 
which mass was swept up (dM/dt). Over time, the dust 
is sputtered and the gas cools. The sputtering of the dust 
occurs on shorter timescales, and thus the integral of the 
IR emission over time will be more heavily weighted to- 
wards the most recently swept up ISM where the small 
grains are still present. 

Consistency between the IR and X-ray analyses can be 
strengthened by adopting the X-ray derived T e and n e t 
(and thus Aa) as fixed parameters of the dust models. 
The work by Hwang et al. (2005) includes a detailed ex- 
amination of a slice across the BEK passing through the 
Compact Knot and the Bar-N regions. Their analysis 
reveals that T e rises by a factor of 2 from the east to 
west side of the Compact Knot and then remains rela- 
tively constant across the Bar. The ionization timescale 
n e t is low (~ 0.5 x 10 11 s cm" 3 ) across the east side and 
center of the Compact Knot, with a sharp jump to higher 
values (~ 2.5 x 10 11 s cm” 3 ) occurring on the western 
edge of the Compact Knot. This jump occurs within 
the region averaged for the IR spectrum of the Compact 
Knot. The ionization timescale remains high across the 
Bar. We have examined these models to determine the 
mean T c and n e t over the three regions corresponding 
to our IR spectra (Table 2). Based on these results, if 
we restrict our grid of IR models to T e = 5.6 x 10® K 
and Aa = 0.0040 - 0.0100 pm (the range allows for the 
jump in n e t) for the Compact Knot, and T e = 5.6 x 10° 
K and A a = 0.0100 pm for the Bar regions, then the 
remaining free parameter of the models, n e , is found 


to increase relative to the best-fitting models to values 
of n e = 5.7 and 8 for the Compact Knot and Bar re- 
gions respectively. However the quality of the fits suf- 
fers. The increased amount of sputtering further reduces 
the contribution of small dust grains, which are respon- 
sible for the short wavelength emission. To compensate 
for the weaker short wavelength emission, the derived 
gas density is higher, which raises the temperature of 
the remaining larger grains. Yet the resulting spectra 
still fall too steeply at short wavelengths (A < 20 pm), 
and the spectra fall too steeply at the longer wavelengths 
(25 < A < 70 pm) as well. When thus constrained, the 
model spectra become too sharply peaked because of the 
lack of the smaller grains with their widely fluctuating 
temperatures. 

4.3. Other Sources 

The IRAC and MIPS images incidentally reveal several 
very distinctive objects which are unassociated with Pup- 
pis A, e.g: (1) A symmetric 16"-radius edge-on galaxy 
at (a, <5) = (8 h 24 m 13 s , — 42°56 / 18") is clearly seen at all 
4 IRAC wavelengths. This galaxy is faintly visible in 
2MASS images at K, H , and (barely) J bands. At 
24 pm, any emission from the galaxy is lost behind the 
brighter emission of Puppis A. (2) The MIPS 24 pm im- 
age contains several compact clusters of point sources 
in the northern half of the field. The brightest ob- 
jects in the clusters are usually detected at 70 pm as 
well. The most prominent of these clusters contains more 
than ten 24 pm sources within a 2' radius, including 
IRAS 008222-4214 and IRAS 08222-4212. Another clus- 
ter (identified by Dutra et al. 2004) and its associated 
reflection nebula produce the extended and extremely 
bright source at (a, S) = (8 h 26 m 26 s , -42°32'40"). (3) 
The bright linear 160 pm cloud, about 15' long at 
(a, <5) = (8 h 22.2 m , — 43°21"), is a dark cloud at visible 
wavelengths. The facts that it is a visible dark cloud 
(Feitzinger & Stiiwe 1984), and that the measured ve- 
locity of its CO emission is +6 km s" 1 (Ramesh 1994), 
compared to ~ +16 km s” 1 for clouds associated with 
the SNR, both indicate that this cloud lies well in the 
foreground of Puppis A. 

5. SUMMARY 

Spitzer MIPS images at 24 and 70 pm show that the 
IR to X-ray correlation within the Puppis A SNR is very 
strong on all spatial scales. The total IR luminosity of the 
SNR is ~ 1.4 x 10 4 L© produced by 0.25 Al© of dust.. The 
ratio of IR to X-ray emission is 6.9, indicating that the 
IR emission is the dominant radiative energy loss of the 
SNR. At 160 pm the SNR becomes lost in the confusion 
of the Galactic plane indicating that there is not a large 
column density of colder dust associated with the SNR. 
At the BEK, IRAC images show weak emission that has 
colors similar to the PAH emission of the general ISM. 
This may indicate the zone where PAHs are destroyed 
in the shock front. Some 5.8 and 8 pm emission asso- 
ciated with optical emission of singly -ionized species is 
also noted in a region south of the BEK. Across the SNR. 
there are only a small number of locations where optical 
filaments have IR counterparts in the broadband MIPS 
images. The most prominent of these is the oxygen rich 
Omega filament. 
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The IRS SL and LL spectra of the East Cloud out- 
side the SNR reveal a typical ISM spectrum, dominated 
by PAH emission features. The longer wavelength ther- 
mal emission of the larger graphite and silicate grains 
suggests an ISRF that is ~ 10 times stronger than the 
local ISRF in the solar vicinity. The spectral map of the 
BEK region does not show PAH features, but does show 
the characteristic silicate band at 20 gm. The external 
grain size distribution must be subjected to sputtering 
to reduce grain sizes by A a « 0.002 pm in order to fit 
the observed BEK spectrum. This implies an exposure 
timescale of ~ 80 yr, which is shorter than timescales 
estimated from prior analysis of the X-ray spectra and 
morphology, indicating a bias towards emission from the 
more recently swept-up dust in a process of extended du- 
ration. In the region of the BEK, the 500 km s -1 shocks 
have destroyed ~ 25% of the mass of the swept-up inter- 


stellar dust. The amount of dust destroyed will increase 
until the post -shock gas cools and sputtering becomes 
ineffective. 
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located at NASA Goddard Space Flight Center. This 
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ditional funding from the NSF through grant AST- 
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